#ifndef DENSITY_H_
#define DENSITY_H_
/*
 * functions.h
 *
 *  Created on: Feb 21, 2013
 *      Author: drewlewis
 */


#include <tiled_array.h>
#include <Eigen/Dense>
#include "Global.h"


void Density_Update(ARRAY2 &R, const ARRAY2 &S, const ARRAY2 &F, std::size_t iter,
                   std::size_t rotation_order = 2);

Matrix
intial_density_guess(const Matrix &S, const Matrix &H, const int occ);

void
Purify(ARRAY2 &R, const ARRAY2 &S);


#endif /* DENSITY_H_ */
